Radiation of caustic beams from a collapsing bullet 
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Collapse of an intense (2+l)-dimensional wave packet in a medium with cubic nonlinearity and 
a two-dimensional dispersion of an order higher than parabolic is studied both theoretically and 
experimentally. The carrier waves are microwave backward volume spin waves which propagate in a 
stripe made from a thin ferrimagnetic film and the packet is a spin- wave bullet. We show that before 
being self- destroyed the bullet irradiates untrapped dispersive waves, which is in agreement with a 
previous theoretical prediction. Since, in addition, the ferromagnetic medium is characterized by 
an induced uniaxial anisotropy, this radiation takes the form of narrow beams of continuous waves 
at very specific angles to its propagation direction. Based on our theoretical calculations we find 
that these beams are caustic beams and the angles are the characteristic spin-wave caustic angles 
modified by the motion of the source. 

PACS numbers: 75.30.Ds, 05.45.Yv, 42. 25. Fx 
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Solitons while perturbed are able to emit small- 
amplitude radiation with frequency detuned far from the 
soliton frequency [l|-[3|. This ability, generally recog- 
nized as striking evidence of soliton's wave nature has 
gained dramatically in importance in recent years with 
the discovery of supercontinuum generation in photonic 
crystal fibres whose applications for frequency comb gen- 
eration in metrology, spectroscopy and imaging are more 
than just impressive (see e.g. [5]). It was shown that dis- 
persion radiation from solitons is responsible for a major 
spectral part of the supercontinuum (see e.g. [4, 6]). 

One has to note that small-amplitude radiation from 
solitons is possible, provided the dispersion relation for 
the medium deviates from the parabolic law uoik) = 
cjo + v{k — ko) + D /2(k — /cq)^ built-in to the standard 
(i.e. parabolic) Cubic Nonlinear Schrodinger Equation 
(CNSE) which governs envelope solitons in optical fibers 
[7|, deep-water solitons J8[, and spin- wave envelope soli- 
tons in magnetic films [9|: 
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Here a is the envelope function for the propagating wave 
packet, z is the direction of propagation, and N is the 
nonlinear coefficient (the term iN\a\'^a is called "cubic 
nonlinearity" and gives rise to the equation name). 

The dispersion law is expressed as a series expansion of 
the frequency uj in terms of the wave number k of plane 
linear waves supported by the medium. The sign of the 
second-order dispersion coefficient D = d'^uj{k)/dk'^ is of 
paramount importance for soliton formation: solitons are 
formed provided DN < [10]. 

Karpman [11] theoretically predicted that in the 
same CNSE framework the fourth-order two-dimensional 
(2D) dispersion uj{ky^kz) with a set of coefficients of 
proper signs should lead to wave irradiation from (2+1)- 
dimensional wave packets a(x, t) as well. He considered 



an axially symmetric wave packet in an isotropic medium. 
In his work and in subsequent numerical simulations [HI 
it was shown that a localized (2+l)D wave packet is de- 
stroyed by radiating energy into a non-trapped circular 
wave. 

To the best of our knowledge this radiation has not 
been observed so far. Interestingly, a similar effect of 
conical Cerenkov radiation from optical bullets in an ax- 
ially symmetrical 3D geometry [13] has attracted a lot of 
attention [14]. The conical radiation does not originate 
from a dispersion of a higher-order but from a higher- 
order nonlinearity. Thus, it can be characterized as an 
effect "dual" to one Karpman predicted. 

In this work we experimentally demonstrate the non- 
trapped radiation from a collapsing bullet originating 
from the higher-order dispersion. This phenomenon is 
observed for a spin-wave bullet propagating in a thin 
magnetic film. Furthermore, we observe two extra fea- 
tures which were not included in the original theory 
[ll| . First, (i) we find that in a medium with uniaxial 
anisotropy this radiation takes the form of two narrow- 
aperture beams. Second, (ii) as the medium is charac- 
terized by a dispersion of backward type, the radiation 
is in the backward direction with respect to the bullet 
propagation direction. This is in contrast to Cerenkov 
conical radiation which is forward [l3| . 

Stable (2+l)D localized nonlinear spin-wave excita- 
tions, termed spin- wave bullets, have been previously ob- 
served in thin ferrimagnetic films of yttrium iron-garnet 
(YIG) magnetized along the propagation direction 
18]. The waves propagating along the field are called 
backward volume magnetostatic spin waves. They are 
waves of backward nature, which means that their group 
velocity and wave vector point in opposite directions. 

Two-dimensional pulses are intrinsically unstable and 
undergo nonlinear narrowing leading to collapse as non- 
linearity overcompensates linear broadening due to 2D 
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FIG. 1. Snapshots of spin- wave packets measured at different 
input powers and at 50 ns time intervals. Arrows show the 
bullet propagation direction. It coincides with the direction 
of the bias magnetic field H and the coordinate axis z. 



parabolic dispersion. Weak magnetic losses in a real mag- 
netic film may balance the nonlinear narrowing. This re- 
sults in a quasi-2D spatially localized bell-shaped wave- 
form called a bullet and ensures its stability for some dis- 
tance of propagation [15]. This regime is well described 
by the (2+l)D parabolic CNSE. The stable bullets are 
observed in a certain range of initial powers for the wave 
packets. For larger input powers waveform collapse is 
unavoidable [15]. The advanced stage of the collapse is 
beyond the limits of applicability of the parabolic CNSE, 
and the specific scenario for destruction of the collapsing 
wave packet ( "collapse scenario" ) should be governed by 
processes of higher-order. 

Our experiment is carried out using a longitudinally 
magnetized yttrium iron garnet (YIG) film stripe which 
is 2.5 mm wide {w = 2.5 mm) and 5 pm thick. The mag- 
netizing field H is 1831 Oe. The spin waves are excited 
by a rf magnetic field created with a 25 pm wide mi- 
crostrip antenna placed across the stripe and driven by 
20 ns long microwave current pulses at a carrier frequency 
of 7.125 GHz. The spatio-temporal behavior of the trav- 
eling spin- wave packets is probed by means of space- and 
time-resolved Brillouin light scattering spectroscopy [19]. 
For relatively small input powers we observe formation 
of quasi-stable (2+l)D wave packets - guided spin-wave 
bullets (Fig. 1, upper row) - reported previously for the 
same geometry [18]. If we increase the power beyond 
the range of bullet stability the wave packet collapses 
(Fig. 1, lower row). The most prominent feature of the 
collapse is the pair of rays irradiated from the packet in 
the backward direction (second panel in the lower row of 
Fig. 1). The left panel in Fig. 2 gives a detailed view of 
the beam emission. One clearly sees that the rays have 
narrow apertures and are directed at weh defined angles 
to the longitudinal axis of the ferromagnetic stripe: the 



FIG. 2. Irradiation of caustic beams from the collapsing bul- 
let. Left panel: experiment. Right panel: simulation. 



value of the angle between the rays is 64°. 

In order to reveal the origins of the irradiation we carry 
out numerical simulations. An original reciprocal-space 
approach is used which effectively accounts for all disper- 
sion orders of the real spin-wave dispersion [ii,|2o|. In 
the reciprocal space CNSE takes the form of a system of 
coupled equations for the amplitudes of the wave packet's 
spatial Fourier components: 
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where a = a{y^z^t) is the spin- wave scalar amplitude and 
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denotes a 2D Fourier transform which is continuous along 
the stripe {kz takes continuous values in the model) and 
discrete in the transverse direction {ky = nir/w^n = 
1,2,3,...). This term is responsible for nonlinear cou- 
pling of Fourier components with different values of n and 
kz. r denotes the coefficient of spin- wave linear damping 
for the medium and fn,k{t) is the Fourier transform of the 
time envelope of the driving microwave field which excites 
the spin- wave packet at the entrance into the waveguide 
medium. Similar to the experiment, fn,k{t) is taken in 
the form of a 20 ns-long rectangular pulse. The carrier 
frequency of this driving pulse is ujo/27t (see Fig. 3). 

The term ujn,k = ^{ky^ kz) is the 2D spectrum of spin 
waves, shown in Fig. 3 as /C;^ -dependence of frequency 
for a family of different discrete values of ky which we 
take here equal to nn/w. The calculated dispersion cor- 
responds to the conditions of our experiment. Physically 
this (n, /c2)-representation describes a family of guided 
modes for a film waveguide. The transverse (?/-)profiles 
of the width modes are given by sm{n7ry/w) with n in- 
dicating each particular mode. One sees that the disper- 
sion slope dujn,k/d\kz \ is negative for all the modes. This 
reflects the fact that the modes are of backward type. 
Consequently, a localized wave packet which propagates 
in the positive direction of the axis z has a carrier wave 
number kz < 0. 

As previously shown for both ID [2] and isotropic 2D 
media [ll|, to ensure radiation from nonlinear wave- 
forms it is important to have a dispersion law for a 
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medium which contains terms higher than parabohc. 
More precisely, the curvature d'^uj/dk'^ should change 
sign along the dispersion curve. From Fig. 3 one sees 
that on one hand for each particular width mode the 
sign of d'^uon.k/ dkl changes from negative to positive in 
the vicinity of /c^ = 0. On the other hand, the spec- 
trum in Fig. 3 is quite dense. Considering this spectrum 
as quasi-continuous one retrieves the 2D continuous dis- 
persion of plane spin waves in ferromagnetic films [21 1 



which is characterized by the variation of dispersion sign 
d^ijjjdk^ as a function of the propagation angle with re- 
spect to the direction of the applied field H (k? = ky-\-kl^ 
(j) = dj:cid.ii{ky/kz)). 

Our simulations confirm the validity of Karpman's idea 
that the higher-order dispersion alone may give rise to a 
non-trapped radiation from a collapsing bullet. The re- 
sult is shown in Fig. 2. (The same parameters as for 
Fig. 3 and standard for YIG films have been used in this 
calculation.) One sees good agreement with the experi- 
ment. In particular, we obtain the same radiation angle 
of 64° as in the experiment. 

The remaining part of the paper is aimed to under- 
standing the extra features (i) and (ii) (see above) of the 
observed radiation and to clarifying the specific radiation 
angle shown in Fig. 2. 

Waves excited by an immobile point source in the 
medium with the dispersion shown in Fig. 3 are prone 
to forming caustics [22|. It is appropriate here to give 
a short introduction to the phenomenon of caustics. In 
an anisotropic medium with a dispersion law co'(k), the 
direction of the wave group velocity = 9a;(k)/9k in- 
dicating the direction of energy flow does not, in general, 
coincide with the direction of the wave vector k. When 
the medium anisotropy is sufficiently strong, the direc- 
tion of the group velocity of the waves in the vicinity of 
a certain wave vector kc may become practically inde- 
pendent of their wave vectors. In this case the energy 
of the waves is channeled along this specific direction, 
which is called caustic direction^ and forms a so-called 
caustic beam. In our recent work [22] we demonstrated 
this phenomenon for spin waves in a single crystal YIG 
film using a quasi-point excitation source whose diameter 
d was smaller than 27r//Cc. In this medium magnetocrys- 
talline anisotropy is negligible; application of an external 
static field H in the film plane imposes a uniaxial sym- 
metry necessary for the formation of caustics. Two pairs 
of spin- wave caustic beams are directed at specific angles 
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One notices that a collapsing spin-wave bullet in a 
magnetically saturated ferromagnetic film meets the size 
criterion d < 27r /kc for a quasi-point source of linear exci- 
tation of caustic waves. Such a source located on the axis 
of the waveguide will excite all waveguide modes which 
have an anti-node on the axis. These are the symmetric 
modes n = 1, 3, 5, . . . For an immobile linear source, the 
frequency of the excited waves is equal to the frequency 
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FIG. 3. Spectrum of the guided modes for the longitudinally 
magnetized film waveguide. The calculation assumes pinned 
spins at the stripe edges. The family of symmetric modes n — 
1,3, ...65 is shown. The mode which is the lowest in frequency 
is the fundamental width mode for the waveguide n = 1. 
Inset: a section of the spectrum close to the frequency cjs/Stt. 
The horizontal blue dashed lines in the main field of the figure 
and in the inset show the carrier frequency cjo/Stt of the input 
microwave pulse. The oblique red lines show the Doppler- 
shifted frequencies of an excitation source which moves in the 
positive direction. The bold sections of these lines show the 
modes responsible for formation of the respective caustics. 



of the source. This condition is indicated by the dashed 
horizontal line in Fig. 3: all the symmetric modes whose 
dispersion lines cross the dashed line will be excited with 
the frequency given by the ordinate of the point of the 
cross-section. It is a short exercise [13] to find which of 
the modes are responsible for caustic formation: those 
modes are shown by the bold part of the dashed horizon- 
tal line. For the conditions of our current experiment the 
angle between the beams in each caustic pair should 
amount to 84° as our calculation based on the theory 
in Ref. [24] shows. Note that this value is significantly 
different from the angle seen in Fig. 2. 

This disagreement is caused by the motion of the bullet 
with the velocity v = 2.5 cm/ps. The moving source is 
described by the red oblique line = cjo + vkz which 
can be considered as the time-space Fourier transform 
of the moving point source (see e.g. Eq. (14) in [23). 
The term vkz is obviously the Doppler frequency shift 
for the excited modes. One sees that now each mode is 
excited with its own frequency. Furthermore, one sees 
that crossings are possible only for positive kz values. 
Importantly, our calculations show that, similar to the 
previous case of the immobile source, the modes which 
satisfy the Doppler shift condition ujn,k{kz) = cjq + vkz 
are also able to form caustic beams. The family of modes 
which are found to be responsible for formation of the 
modified caustics is shown by the bold section of the red 
oblique line. 

The modified caustic angles are obtained by consider- 
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FIG. 4. Slowness curves for the immobile and the mobile 
excitation sources (a) and the respective angles (b). In (a): 
bold solid lines are the slowness curves for the frequency given 
by the horizontal dashed line in Fig. 3. Thin solid lines are 
the slowness curves for the Doppler-shifted frequencies. The 
red bold dashed line connects zero- curvature points of the 
thin solid lines. Dark blue arrows are normal to the thick 
solid line and show the caustic directions for excitation by an 
immobile source. Grey arrows are normal to the family of the 
thin solid lines and show the energy flow directions (DEF) for 
the moving source. In (b): the arrow indicates the direction of 
source motion. Thin blue solid lines: caustic directions for the 
immobile source. Red dotted lines: two DEF for the moving 
source. Green dashed lines: directions along the observed 
caustic beams. 



ing the family of slowness curves for the Doppler-shifted 
frequencies. A slowness curve is a constant-frequency line 
in the {ky^ kz) plane (Fig. 4a) calculated for a film which 
is continuous in both in-plane directions . The group 
velocities of all plane waves which exist at the respective 
frequency are directed perpendicular to this curve. For 
this reason the direction of the energy flow (DEF) [24] is 
given by the normal to the slowness curve at the point 
where its curvature is zero. A range of frequencies for 
which the slowness directions are close to each other will 
contribute to formation of the modified caustics. The 
thin solid lines in Fig. 4 are the slowness directions for 
the frequencies which correspond to the crossing points 
of the oblique red line with the dispersion lines in Fig. 3. 
The bold dashed line in Fig. 4 runs across zeros of curva- 
ture for these curves. One sees that normals to the thin 
solid lines are collinear. This suggests that the normals 
define the DEF directions. The rigorous values for the 
angles between these directions and the bullet velocity 
are ±136.5° (Fig. 4b). 

First, one sees that the angles exceed 90° which reflects 
the fact that the carrier waves are the waves of backward 
nature. (Note that the Doppler frequency shift for back- 
ward waves is anomalous [28, 29].) Second, one notices 
that the angle (j)M = 87° between the two DEF is in- 
creased by just 3° with respect to <pi = 84°. 

Furthermore, for a moving source the directions along 
the formed beams do not coincide with the directions 
of the energy flow. It is instructive to define DEF 
via the angle of group velocity v^. One finds (j)M = 



2 diTctdin{\vgy/vgz\). The direction along a beam is related 
to the direction of by a Galilean transformation. As 
a result the angle between the two backwards-irradiated 
beams is given by 2 arctan( | Vgy \ / \ Vgz —v\). Using this for- 
mula one obtains a value of 64° which is in the excellent 
agreement with the results shown in Fig. 2. 

In conclusion, we studied the collapse scenario for an 
intense two-dimensional wave packet propagating in a 
ferromagnetic medium and experimentally showed that 
before being self-destroyed the wave packet irradiates un- 
trapped continuous waves. Our simulation proves that 
the observed effect is the fundamental phenomenon of 
wave emission from collapsing packets in materials with 
cubic nonlinearity and 2D-dispersion of an order higher 
than parabolic. In addition, since the ferromagnetic-film 
medium used in the present study is also characterized 
by a uniaxial anisotropy, the observed radiation takes 
the form of narrow-aperture beams of continuous waves 
at very specific angles to bullet's propagation direction. 
These angles are larger than 90 degree, since the under- 
lying wave dispersion is of backward type. This specific 
effect of radiation of untrapped waves along the caustic 
directions may exist for other 2D-media media as well, 
provided uniaxial anisotropy of dispersion is available or 
induced in the medium. 
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